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INTRODUCTION 



Recently, iKashlinskv et all (l2008l. 



lAtrio-Barandela et al. ( 2010l ): IKashlinskv et all 



analyzed CMB fluctuations on directions of X-ray galaxy 
clusters and found a large bulk flow with speed wdf ~ 10^ 
km/s to zoF ~ 0.3, toward direction {lo,bo) ~ (290°, 30°). 
This extraordinarily large bulk flow (the so called dark 
flow) is a severe challenge to the standard ACDM paradigm 
and has fundamental implications on the topology of 
the universe, t he inflationary scenario and the nature 
of gravity ( e. g. ICarroU et al. 2008*: lAfshordi et al.l 



Chang et al. 



i.gJUM 
l200g : 



Khourv fc WvmanI l2009l : iKashhnskv et all |2010| and 



ABSTR ACT 

Recently IKashlinskv et al.l (j2008l l2010l) reported a discovery of a ~ 10"^ km/s bulk 
flow of the universe out to z ~ 0.3, through the dark flow induced CMB dipole in 
directions of clusters. We point out that, if this dark flow exists, it will also induce 
observable CMB temperature fluctuations at multipole £ ~ 10'^, through modulation of 
the inhomogeneous electron distribution on the uniform dark flow. The induced small 
scale kinetic Sunyaev Zel'dovich (SZ) effect will reach ^ at multipole 10'^ ^ ^ ^ 

10'', only a factor of ~ 2 smaller than the conventional kinetic SZ effect. Furthermore, 
it will be correlated with the large scale structure (LSS) and its correlation with 
2MASS galaxy distribution reaches 0.3/zK at £ = 10^, under a directional dependent 
optimal weighting scheme. We estimate that, WMAP plus 2MASS should already be 
able to detect this dark flow induced small scale kinetic SZ effect with ^ 6(t confldcnce. 
Deeper galaxy surveys such as SDSS can further improve the measurement. Planck 
plus existing galaxy surveys can reach > 14cr detection. Existing CMB-LSS cross 
correlation measurements shall be reanalyzed to test the existence of the dark flow 
and, if it exists, shall be used to eliminate possible bias on the integrated Sachs- Wolfe 
effect measurement through the CMB-LSS cross correlation. 

Key words: (cosmology:) large-scale structure of Universe: cosmic microwave back- 
ground: theory: observations 



ture fluctuations. If the amplitude of the dark flow reaches 
the reported value of ~ 10^ km/s, it shall be detected by 
combining existing CMB and galaxy measurements. 

If the matter distribution where the dark flow resides 
is homogeneous, the dark flow only induces a CMB dipole. 
However, from galaxy surveys, we know that the local uni- 
verse is_strongl^_^nhomogeneous in the density distribution 
(e.g. iTegmark et"ai] l2004l l. Given these inhomogeneities, 
even an uniform dark flow can induce small angular scale 
temperature fluctuations in the CMB sky, through exactly 
the same mechanism of the inverse Compton scattering to 
genera te the conventional kinetic Suny aev Zel'dovich (SZ) 
effect l|Sunvaev fc Zeldovichll 19721 . Il980l ) , 



Mersini-Houghton fc Holman 



200£ 



200£ 



references therein). 

The above dark flow measurements ar e under scrutiny 
(|Keisleill2009l : lAtrio-Barandela et al.|[201ol ). It also requires 
indepen dent confl rmations, from dir ect (iFe ldman et al.l 
[2009; Watkins et al. 2009) and indirect (|Lavau"x ct al. 2oTdi ] 
velocit y measurements of n earby gala:xies and nearby super- 
novae ll Gordon et al ] |2008l ). In this paper, we point out that, 
besides the dipole from which the dark flow is inferred, the 
dark flow also induces small angular scale CMB tempera- 
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Here, cos 6 = tt-df ■ n is the cosine between the dark flow 
direction jidf and the direction h. zdf is the edge of the 
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dark flow, — he{l + Se) is the 3D free electron number 
density, fie is the mean number density and Se is the over- 
density. Throughout the paper we adopt the fiducial value 
f2i,/i — 0.031 and thus neglect the prefactor f2b/!,/0.031 in 
Eq.Hl X = x/i^/Ho) is the dimensionless comoving distance 
in unit of the Hubble radius c/Hq. Throughout the paper 
we adopt a flat ACDM cosmology with ~ 0.27, f^A = 
l-nm,ilb = 0.044, as = 0.84 and h = 0.71. 

The flrst term in the last expression of Eq.[T]is a dipole 
term. However, the second term is not, despite its dipole-like 
prefactor cos^. The electron over-density 5e{n) has a com- 
plicated directional dependence and is clustered at scales 
< 100 Mpc. It is this density modulation from all free 
electrons producing the small scale kinetic SZ effect, the one 
that we point out and investigate in this paper. It clearly dif- 
fers from the CMB dipole induced by the Earth motion and 
the CMB dipole induced by the dark flow of uniform elec- 
tron distribution, both are lacking of the &e modulation. It 
also differs from the conventional kinetic SZ effe ct, which 
is fur ther modulated by the non-uniform velocity l|Vishniacl 
ligsTl ). Later we will find that the two kinetic SZ effects have 
different clustering behaviors for this reason. 

This dark fiow induced small scale kinetic SZ effect 
leaves unique imprints on the CMB sky. We will estimate 
its auto power spectrum and its correlation with the large 
scale structu re (LSS). The exa ct value of hdf is highly un- 
certain (Kashlinskv et al]l2010l ). Throughout the paper, we 
will adopt a fiducial value wdf ~ 10^ km/s. The results 
presented in this paper can be scaled to other value of vdf 
straightforwardly. Given -udf = 10^ km/s, we find that it 
is promising to extract the dark flow component through 
CMB-LSS correlations in an unbiased manner, if the CMB 
maps are properly weighted with a directional dependent 
weighting factor. We estimate that existing data may al- 
ready allow for > 6(T detection of this effect and provide 
independent test on the dark flow scenario. 



2 THE AUTO POWER SPECTRUM OF THE 
DARK FLOW INDUCED KINETIC SZ 
EFFECT 

The dark flow induced kinetic SZ temperature fluctuation is 
correlated at small scales due to clustering of the underlying 
electron overdensity 5e- We derive the resulting auto power 
spectrum averaged over the survey area (the appendix ^K\ . 
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The expression in the integral has adopted the well known 
Limber approximation. Ae(fc, z) is the electron number den- 
sity power spectrum (variance) at scale k and redshift z. 
Throughout the paper we approximate it as = Aj^. 
Here A^ is the matter power spectrum (variance). Since 
the measured matter distribution i n the nearby univers e 
agrees with the standard ACDM (|Tegmark et all l2004h . 
it allows us to adopt the CMBFAST transfer function 
l|Seliak fc Zaldarriaga|[l99^ . The nonlin ear power spectrum 
A^ is calculated by the halofit formula jSmith et al.ll200"3 ). 



The suppression factor fr is 



(3) 



The average is over all directions in the survey area. Clearly 
/t is a direction dependent quantity and varies from sur- 
vey to survey. Namely, the dark flow induced kinetic SZ ef- 
fect is statistically anisotropic, unlike other isotropic CMB 
components. Eq.[3]is derived under the small angle approx- 
imation. Its accuracy is of the order {Ad)^/4 ~ (tt/^)^ <C 1 
when £ ^ 10, where A9 is the typical angular scale involved, 
A^ ~ 2n/£. Refer to the appendix ^X]for more details. For 
a full sky survey /t = 1/3. 

This characteristic directional dependence can in prin- 
ciple be applied to separate the dark flow induced kinetic SZ 
effect from isotropic components (primary CMB, the ther- 
mal SZ effect, the conventional kinetic SZ effect, etc.) in the 
auto correlation measurement. Alternatively one can weigh 
the CMB temperature ffuctuations by a directional depen- 
dent function W{n). It can be desired to amplify the dark 
flow component with respect to other components. In this 
case. 



/t = {cos'' ew^ih))s 



(4) 



Cdf(^) shows a number of difference in the £ depen- 
dence comparing to that of the conventional kinetic SZ ef- 
fect (Fig. [l}, the major reason is that the dark flow in- 
duced kinetic SZ fluctuation is linear in density fluctua- 
tions while the conventional one is quadratic through the 
interplay between density and velocity inhomogeneity. An- 
other reason is that the dark flow induced kinetic SZ effect 
comes from much lower redshift than the conventional ki- 
netic SZ effect. Since the typical redshift of the dark flow 
induced kinetic SZ effect is proportional to zof, the peak 
moves to smaller scale (larger £) when zof increases. For 
the fiducial value of wdf = 10^ km/s out to z-df = 0.3, 
the dark ffow induced kSZ power spectrum, averaged over 
the whole sky with W = 1, peaks at ^ ~ 2000 with an 
amplitude LQ/jK^, a factor o f ^ 2 smaller than the conven- 
tional kinetic SZ effect fe.g. [Ma &. FrvlpOo3 : IZhang et all 



l2004lHernandez-Monteagudo fc Hol2009l ') at the same scale 
Overall, this dark flow induced kinetic SZ effect is non- 
negligible, ranging from ~ 100% of the conventional one 
at ^ = 200 to ~ 30% at ^ = 10^. 

The SPT collaboration measured a combined SZ power 
spectrum (tSZ-f0.46 x kSZ) at 150 GH z band and £ = 3000 
to be 4.2 ± 1.5/iK^ (|Lueker et al.ll2009l ). This measurement 
is a factor of ~ 2 lo wer than the predic t ed thermal SZ effec t 
with (78 = 0.8 fe.g. IZhang et al.fl2002l : iLueker et al.ll2009l ). 
but consistent with rece nt WMAP and ACT mea surements 
(Komatsu ct al. 2010; The ACT Collaboration e t al. 201(i). 
The dark ffow induced kinetic SZ effect is ~ 0.3/^K^ at the 
analyzed SZ sky where fr — 0.14. Given its sub-dominance 
to the conventional kinetic SZ effect and the thermal SZ ef- 
fect (Fig. [T|, it is difficult to perform a robust test of the 
reported dark flow against existing SZ power spectrum mea- 
surements. Future large-area SZ measurements at 217 GHz, 
the null of the thermal SZ effect, may be able to infer its 
existence through the directional dependence in the fr pref- 
actor. 
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Figure 1 . The auto correlation power spectrum of the dark flow 
induced kinetic SZ effect. Th e primary CMB data points are the 
seven years WMAP result llKomatsu et alj [2OIO; Larson et al. 
I2OICI) . The SZ measurement at £ = 3000 is from Lueker et al. 
1I2OO9I '). The thermal SZ result is the one in jZhang et al.j 12002), 
scaled from the original erg = 1.0 to eg = 0.8 with the (t| scal- 
ing and scaled from the Raleigh-Jeans regime to 150 GHz. The 
kinetic SZ result is the hom ogeneous kine t ic SZ power spectrum 
predicted from the model of IZhang et al.l l|2004l ) . The four solid 
lines are for the dark flow models with z^p = 0.1,0.2,0.3,0.4 
from bottom up. We have adopted v-Q-p = 10^ km/s. The shown 
power spectra is for a full sky survey, with th e suppression factor 
/y = 1/3. The SPT sky l|Lueker et al.l l2009h has /t ~ 0.14, so 
the dark flow contribution to the measured SZ effect at £ = 3000 
is ~ 0.3^K^, sub-dominant to other components. 



3 CROSS CORRELATION WITH GALAXIES 

The dark flow induced kinetic SZ effect is also correlated 
with tracers of the large scale structure at 2; < zuf- The re- 
sulting cross power spectrum between the dark flow induced 
kinetic SZ effect and the galaxy distribution is 

«2/- 
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lO^km/s 



(5) 
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fig is the galaxy distribution function normalized that 
J ng{z)dz = 1. We approximate the electron-galaxy cross 
power spectrum A^g = &gA^ where bg is the galaxy bias. 
The suppression factor 



fTg = {cos eW in)) ^ 



(6) 



Here the average is over the overlapping sky of CMB and 
galaxy surveys, frg is also direction dependent and varies 
from survey to survey. Especially, fxg ~ with W — 1 
for fuU sky surveys. Thus for WMAP-f2MASS, we have to 
choose W 7^ 1 to avoid the cancellation. A natural choice is 
W{n) — cos 9, for which fxg = 1/3. A nice property about 



Figure 2. The cross correlation power spectrum between the 
dark flow induced kinetic SZ effect and a full sky galaxy survey 
of 2MASS-like. A directional dependent weighting W = cos 9 is 
adopted. Since most 2MASS galaxies reside at 2 < 0.1, the cross 
correlation varies little for bulk flows beyond 2 > 0.1. For this 
reason, the three curves of 2df = 0.2, 0.3, 0.4 are barely distin- 
guishable. 



this choice is that {W)s ~ and thus isotropic components 
in CMB maps (primary CMB, the thermal SZ, conventional 
kinetic SZ effect, dusty star forming galaxies, etc.) do not 
bias the cross correlation Q 

The cross correlation signal between the dark ffow in- 
duced kinetic SZ effect and 2MASS galaxy distribution is 
shown in Fig. [ 2l Th e galaxy distribution is adopted from 
lAfshordi et al.l ( 20041 ) . along with the galaxy bias bg — 1.18. 
Since the galaxy distribution peaks at z < 0.1, the peak of 
the cross power spectrum is at lower £ than that in the auto 



power spectrum. For a dark flow with vof = 10 km/s out to 
zdf = 0.3, the cross power spectrum peaks at ^ ~ 1.3 x 10^ 
with amplitude 0.32/iK. 

Given the existence of this cross correlation, it is inter- 
esting to ask whether it will impact the interpretation of 
existing CMB-LSS cross correlation measurements, whose 
primary goal is to detect the inte grated Sachs- Wolfe effect 
(ISW) (|Giannantonio et al.l (|2008l ) and references therein). 
Since fxg = with W = 1 for full sky surveys, this 
dark flow induced kinetic SZ effect is unlikely to signifl- 
cantly bias the integrated Sachs- Wolfe (ISW) measurement 
through WMAP-f2MASS (NVSS), unless the CMB masks 
causes significant deviation of frg 7^ 0. For surveys with 
partial sky coverage, frg 7^ in general. So the ISW mea- 



^ Realistic correlation analysis masked the galactic plane, so 
/rg ^ even for full sky surveys like WMAP-I-2MASS. This is- 
sue is important for actual data analysis. But for the theoretical 
study presented in this paper it is safe to neglect it. 
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Figure 3. The predicted S/N of WMAP-2MASS 
(Planck+2MASS) cross correlation designed for the dark 
flow detection, for which an unbiased optimal directional de- 
pendent weighing scheme is applied. WMAP+2MASS is able to 
detect the bulk flow at ~ 6cr confidence level. Planck-|-2MASS is 
able to improve the S/N by a factor of ~ 2. Since most 2MASS 
galaxies reside below z = O.f, it is inefficient to detect the bulk 
flow beyond z = 0.1, explaining the plateau in the S/N curves. 
Cross correlating WMAP and Planck with deeper galaxy surveys 
such as SDSS, LAMOST and BOSS will improve the detection 
at z > O.f. Thus existing surveys (WMAP, 2MASS and SDSS 
esp ecially) are alrea d y abl e to d etect the dark flow reported 
by iKashlinskv etID (2OO8I . I20f0h . Since the S/N oc vuf, the 
cross correlation measurement of existing surveys can detect 
dark flow with speed ~ 500 km/s at > 3cr confidence level and 
Planck-|-2MASS can improve the sensitivity to ~ 200 km/s. 



Here, C'^^^ is the power spectrum of primary CMB and 
C^^^'N is the one of other components including the in- 
strumental noise. Since Cdf < c'^^^ _|_ c"<sz ^ 

safely neglect the contribution of Cdf in (7Cmb,n ^j^^ 
galaxy power spectrum and Cg,jv is the associated measure- 
ment noise, Cg,]v = Anf^j^y/Ng where Ng is the total number 
of galaxies, /sky is the fractional sky coverage. 

We not only want the weighting to be unbiased ({W)s = 
0), but also want it to be optimal such that the measure- 
ment error is minimized. It turns out the unbiased optimal 
weighting is the solution to a 2nd integral equation derived 
in the appendix ^B] 

The forecast for WMAP (PLANCK)-2M ASS cross cor- 
relation measurement is shown in Fig. |3] We take the im- 
pact of galactic mask into account and adopt /sky = 0.7 
to evaluate the cosmic variance and shot noise. The unbi- 
ased optimal weighting W for this survey configuration is 
not trivial to derive. However, for the error estimation pre- 
sented here, it is safe to adopt the one corresponding to 
the case of full sky coverage, for which we have the analyt- 
ical solution W{n) = cos 9. Under such weighting we have 
{W^)s = frg = 1/3. We find that WMAP plus 2MASS is 
able to measure the dark fiow, if it extends to z ^ 0.1 with 
an amplitude 10'' km/s, at ~ Qa. Due to the resolution of 
WMAP, it will miss the peak of the correlation at ^ ~ 10''. 
On the other hand, Planck, with better angular resolution, 
will be able to well capture the peak correlation and is thus 
able to improve the measurement to ~ 14a. 

Since most 2MASS galaxies locate at z < 0.1, 2MASS 
is inefficient to probe the dark flow at z > 0.1. This ex- 
plains the S/N plateau at zdf > 0.1 (Fig. |3]). Deeper sur- 
veys such as SDSS are more suitable for this purpose and 
can improve the overall S/N significantly. We expect that 
combining WMAP and existing galaxy surveys are already 
able to detect the dark fiow ai z < 0.3 with the claimed am- 
plitude and depth at > 6a level. Planck will further improve 
it to > 14a, for vuf = 10'' km/s. Since the S/N oc udf in 
the cross correlation measurement, this implies that, Planck 
plus existing galaxy surveys are able to detect dark flow with 
amplitude as low as 200 km/s at > 3a. 



surements based on these surveys are biased by this com- 
ponent. Interestingly, depending on the direction of survey 
area, this dark flow induced kinetic SZ effect may increase 
or decrease the measured correlation strength. Given the 
comparable correlation strength of the two effects, this issue 
shall be taken into account of future data analysis to avoid 
bias in dark energy constraint, especially for low redshift 
galaxy surveys with partial sky coverage like SDSS. 

Thanks to its characteristic directional dependence, at 
least in principle we are able to measure the dark flow 
through CMB-LSS cross correlation in a way unbiased by 
other CMB components such as the ISW effect. To do so, 
we need to design a weighting scheme such that {W)s ~ 
to eliminate isotropic components in CMB maps (primary 
CMB, the thermal SZ, conventional kinetic SZ effect, dusty 
star forming galaxies, the ISW effect, etc.). The statistical 
measurement error under such weighting is 

ACrg _ / l + (ly^)s(CCMB + CCMB,N)(C^ + C3,iv)/C|, 

Crg V 2£A£Uy ^ > 



4 SUMMARY 

We have pointed out the existence of the dark flow induced 
small scale kinetic SZ effect and estimated its amplitude. As 
a potentially non-negligible component of CMB temperature 
fluctuations, it impacts cosmology in at least two ways. 

It enables a useful i ndependent check on existin 



dark flow measurements llKashlinskv et al. I2OO8I. I2OO' 



lAtrio-Barandela et al.l bOlOinKashlinskv et al. I I2OIOI ). The 

direction-weighted CMB-LSS cross correlation measurement 
proposed in this paper should be able to detect a dark 
flow with an amplitude of > 500 km/s at > 3a level for 
WMAP-I-2MASS. Planck-|-2MASS can improve the detec- 
tion threshold to 200 km/s. Given this sensitivity, it will 
allow for a strong test of the existence of the dark flow. 

• A dark flow of the reported amplitude ~ 10"' km/s can 
significantly bias the ISW measurement through the CMB- 
LSS cross correlation, if the survey area is sufficiently close 
to the direction of the dark flow or the opposite of it. Hence 
it can signiflcantly bias the dark energy constraint based on 
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the ISW interpretation. Existing data shall be reinterpreted 
to avoid such bias. 
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APPENDIX A: DERIVING THE SUPPRESSION 
FACTOR 

We first subtract the dipole mode in Eq. [T] and rewrite the 
rest as 



Odf = cos 9 X y{n) ; y{n) oc / 5^(1-1-2) dx 



(Al) 



The expectation value of the angular correlation between 
two fixed directions hi and n2 is 



(9Dp(ni)eDp(n2)> 



coseicos02(y{hi)y{n2)) (A2) 



' + o{'-f) 



{y{ni)y{h2)) 



Here, 9 = (Oi + 62) /2 and 9i2 = 9i - 92. For smaU an- 
gular separation |A^| <^ 1 {hi ■ h2 = cosA^), we have 
I6I12I 5? \A9\ < 1. So the term 012/4 < 1 and can be safely 
neglected. Averaging over all pairs with the same A9 in the 
survey volume and using the fact that {y{hi)y{h2)) should 
be isotropic, we have 

(6DF(ni)6DF(n2))s ^ {cos'^ 0)s{y (hi) y(h2)) . 

Fourier transforming the above equation, we obtain Eq.[2]& 
Eq. Oand recognize fr = (cos^ 9)s for the case of W — 1. 
For a general W, fr = (cos^ 9W^)s can be derived following 
the same procedure. 



APPENDIX B: THE UNBIASED OPTIMAL 
WEIGHTING FUNCTION 

The optimal weighting W(h) shall minimize {Wg)/f^g oc 
{W^)s/{cos9W)l, under the constraint {W)s = 0. By the 
Lagrange multiplier method, this corresponds to minimize 



(cosevy)! 
This requires 



MW}s 



W - 



{cos 9W)s 



[cos 6* - (cos 61) s)] 



(Bl) 



(B2) 



For a general configuration of sky survey area, it is non- 
trivial to solve the above second order integral equation. 
However, for the full sky coverage, since {cos9)s = 0, one 
can easily find the solution to be = cos 9 = huF ■ h. This 
is what we adopt to estimate the expected WMAP-I-2MASS 
cross correlation S/N. 



